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1 Introduction 



In non-relativistic dynamics time has a characterization of its own which 
distinguishes it sharply from the space coordinates of configuration space. 
However, this difference can be simply removed at the formal level by going 
to the parametrized form of dynamics where time is made to depend on a 
parameter r in as much as the coordinates q % do. One is thus led to deal 
with a set (g*(r),t(r)) in which the identification of time versus coordinates 
appears more as a matter of convention than as a matter of significance 
from the point of view of the dynamical system under study. Even though, 
time still keeps a particular role from the physical point of view. Time is 
experienced by the observer as well as by the system. This is more evident 
in the transition to quantum mechanics, where time -as opposed to position- 
can not be viewed as a property of the system under scrutiny. 

There is a way out from this situation as shown in ref [Q, whose authors 
show how to deal with and solve the question at what time? in quantum 
mechanics in one space dimension by introducing a suitable time operator, 
and obtaining the associated time representation. The outcome is the emer- 
gence of a x «-> t equivalence in quantum mechanics in much the same way 
as there is one in classical mechanics. The question at what time? joins the 
question at what position? as answerable not only experimentally, but also 
within the realm of the quantum mechanical formalism. 

In special relativity time is obviously q°, and it seems the question at what 
time? would be addressed in relativistic quantum mechanics in a simple and 
direct way: explicit covariance should rule the presence of q° along with 
the space components q l to form a Minkowski space fourvector g M . There 
should be no telling difference between the time and the space components 
of q, mainly taking into account that -in contrast to the non relativistic 
case- they get entangled by Lorentz transformations. One could be led to 
believe in the existence of a space-time position operator, a four- vector, whose 
components should transform covariantly under the Lorentz group. This 
object should address simultaneously the two questions when? and where? 
seemingly unrelated in the non relativistic case. It is well known that this 
object has never been constructed. In the instant form of dynamics, i.e. 
refering the operators to their values at some instant of time, one can employ 
a three-vector operator -the position operator 0- to answer the question 
where? This operator not only lacks explicit covariance, it also lacks a time 
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component. The cause of these deficiencies can be traced back || to the 
reparametrization invariance of the action of the relativistic particle 



which translates into evolution (along r) generated by a Hamiltonian H = 
p 2 — m 2 = 0. Since the Hamiltonian is constrained to vanish, the r evolution 
is a gauge transformation. In the canonical approach one chooses a solution 
to the constraint, i.e. by putting p° = ^Jp 1 + m 2 , and "fixes the gauge" 
by setting the evolution parameter to be the physical time. A priori there 
is no room left for the question when? as there is no freedom left for a 
time operator differing from the time parameter q . This is a bonus from 
another point of view: demoting q° to the role of a parameter one evades the 
difficulty of a Hamiltonian unbounded from below in the same way as in the 
non-relativistic case. The lack of positivity of the density j° of the solutions 
of the Klein Gordon equation also plays a role here. It brings about particle- 
antiparticle pairs, etc. and the untenability of the one particle interpretation. 
From here on, the true variables are field configurations, to whom q°, along 
with the space coordinates q l , are mere parameters. However, the case of the 
relativistic particle we are analyzing here is of intrinsic interest; it serves to 
set up the basis for the particle interpretation of quantum field theory, and 
also as a guideline to use @] in the construction of the quantum formalism of 
the gravitational field. Analyzing issues of time for the relativistic particle 
may prove valuable in transforming that formalism in a theory or, at least, 
may throw some light on the issues of time in quantum gravity ||. This 
paper focuses on the relativistic particle. In Sect. 2 we summarize the 
results of the canonical formalism, In section 3 we generalize the treatment 
of Ref. Jl| to the free relativistic particle, Sect. 4 contains the generalization 
to three space dimensions and Sect. 5 is devoted to questions of orthogonality 
and completeness. Finally, in section 6 we discuss some issues raised by the 
interpretation of the formalism and some speculations about the applicability 
to quantum gravity. 

2 Canonical formalism 

Here we will focus our attention onto the physical Hilbert space 'Hkg of 
the positive energy solutions ijj(x) for the Klein Gordon equation ||, with 




(1) 
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the understanding that negative energies will be reinterpreted in terms of 
antiparticles. In configuration space where the Klein Gordon equation reads 
(□ + m 2 )ip(x) = 0, the positive energy solutions are of the form: 

ijj{x) = (2tt)- 3/2 / d A ke- ikx 5(k 2 -m 2 )6(k°)V(k) 



with a scalar product: 



( 27r )-3/2 J _^ e -^W*°-^(£) ( 2 ) 



where u(k) — \Jk 2 + m 2 . We will follow the conventions of [0 denoting with 
uppercase letters the wave functions in momentum space, leaving the lower 
case for configuration space functions. 

To answer the question "What is the probability of finding the particle 
at the point x at time x°T with the above scalar product, we need to find 
a Hermitian position operator and find its eigenfunctions ipg tX o. Then, the 
probability amplitude in x for finding a particle at x at time x° = q° is 
i^x,x°i 0), where 0(g) is the wave function giving the state of the particle. As 
shown by Newton and Wigner [[J the position operator is 

3 (4) 

In our notation, k will represent p in momentum space, while Q's and p's 
will denote operators, unless specified otherwise by the word "classically" , in 
which case they will denote classical dynamical variables . The eigenstate of 
the position operator localized at the point x at t = is 



^ j0 (k) = {2-K)^' 2 puj{k)e- tk3 (5) 

In general, given a particle in the state at t — 0, the probability ampli- 
tude to find it at the position x at t = is given by 

(%o, 0) = (2tt)- 3 / 2 / ^f" S y/Mk)Hk) (6) 
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The components of the position operator are in involution and commute 
canonically with the momenta 

[Q\Q j ]=0, [Q\p i }=i5 i i (7) 

under rotations and space translations Q behaves as a three vector. It also 
evolves like the position of a particle should do, namely 

The Heisenberg position operator at time t can be obtained by integrating 
this equation 

Q{t) = Q + -t^t (9) 
uj(p) 

We now would like to invert this equation to get an operator for the time-of- 
arrival of the relativistic particle following the proposal of JTJ. 



3 Time-of-arrival in one space dimension 

The special role played by time has been the source of controversy since the 
early days of quantum mechanics. The search of the various time operators 
and the analysis of the associated time-energy uncertainity relations was 
the subject of a number of works (see the bibliography in ref. ||), whose 
outcome was that quantum mechanics can not accomodate a time-of-arrival 
operator. This has been refuted recently in ref. || where, in addition, an 
average value for this quantity is explicitly obtained for one space dimension 
in terms of the current density of the particle. This is framed in a wealth 
of recent works devoted to the issue of time in quantum mechanics -see 
ref. M and the bibliography contained therein- with special emphasis on 
the tunnelling times, a question of fundamental and practical implications. 
Here, we are interested in the characterization of the time-of-arrival as one 
of the properties of the system under study as in ref. |]], i n other words, 
we need to go one step further and to obtain an associated operator to be 
able to analyze and give an interpretation to this property in the quantum 
formalism. This is necessary for our results to be of value for the quantum 
formalism of the gravitational field where, as said in the introduction, time 
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has to be considered as a property of the system under study. For the sake 
of simplicity and also to connect with the non-relativistic one-dimensional 
case studied in we begin by considering the case of one space dimension. 
Then we can rewrite (^) as 

Q{t) = Q + J f 1 t (io) 

and the time-of- arrival at the position Q(t) = X would be given by a suitable 
ordering of the operator 

Q\X)c{X-Q)^- (11) 

where the simbol ~ is employed to mean equal apart from ordering. Now, 
Q°(X) can be given simply in a form that goes to the operator T(X) of [|TJ 
in the non-relativistic limit: 



The eigenfunctions of this operator 

Q°(X)^ TtX (k)=T^ TtX (k) (13) 

are given by 

* TtX (k) = aVke'MW-W (14) 

where a is a normalization factor. Multiplying by the phase factor exp (—imT) , 
these functions give the eigenfunctions of Q in the non-relativistic limit. We 
will not make distintions between right (k > 0) and left moving (k < 0) 
particles here, as these have a meaning for one space dimension only and 
we want to study the 3-D case, where opposite directions can be connected 
continuously. 



4 Three space dimensions 

A new feature appears in three space dimensions that was not present in 
the case studied above. The space of "detected" states is a subspace of the 
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Hilbert space Hkg of positive energy solutions to the Klein Gordon equation. 
This comes about because in the 3-D case the evolution equations that we 
have to invert to obtain the time-of-arrival is the set of three equations 
depending on a unique parameter t. To be compatible, they have to satisfy 
the constraint 

C=(Q-X)Ap = (15) 

where the "point-of-arrival" X plays the role of a parameter. Classically, 
these constraints mean that the angular momentum of the particle is X A p, 
so that X is a point in the particle trajectory, or simply that the angular 
momentum about X is zero. In quantum mechanics there are obstructions to 
imposing simultaneous values to different components of the angular momen- 
tum. At first sight, the best one can do is to constrain L 2 and a component 
of the angular momentum, say L 3 , to have definite values given from X A p. 
However, this is not the case here, as we are equating the components of the 
angular momentum to an operator X A p, in such a way that the constraints 
form a first class system. Classically, Eq. (fl5|) plays the role of a set of 
first class constraints in the hamiltonian formalism that we have to quantize 
following the method of Dirac. Now, the total hamiltonian is 



F = yy + m 2 + A a C a (16) 

where 

C a = e abc (Q - X) b p c (17) 

and the p's and Q's are the dynamical variables to become operators after 
quantization. It is straightforward to show that 

{C a ,Cb} = e a b c C c , {C a , H} = e a b c \bC c (18) 

Therefore, we have a true first class system, a different one for each vector 
X. 

There seems to be additional difficulties in that the eigenvalues of L 2 and 
L3 are integer numbers while the constraint will assign to them a continuous 



spectrum. Actually, this is not the case |10| because, even if the constraint 
can be written in the form L = X A p, this will not hold as an operator 
equation, nor the states on wich it will be satisfied will be eigenstates of 
neither L« nor X Ap. Now, the detected subspace H K q can be given simply 
as that spanned by the functions ^( x \k) of the form 
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ftgg = = e-^q(k,X)} (19) 

where ty(k,X) represents an arbitrary function of the modulus of k and of 
X. If we now require invariance under translations, we have to drop the 
dependence of ^(k, X) on X. In this case we can say that the Hilbert space 
1~L < kg is obtained from Ti^G by a translation of amount X. 

We are now prepared to study Q°(X), the time-of-arrival at a point X in 
the 3-D space. Classically, it is given by inverting the equation of motion: 

Q°(X) = ^-(Q-X)-p, (20) 
which is a first class dynamical variable {Q°(X),C a } = 0. In the Hilbert space 

(X) 

H K q the operator equation of motion has to be rewritten with t replaced by 
the operator Q°(X) and Q(t) by the detector's position X 

X-Q-JL Q °(X) = (21) 
uj(p) 

It should be an identity, with the operator Q° being such to annihilate the 
left hand side. By vector product of the above equation by p we obtain 
the constraints that are already satisfied in the detected subspace. Scalar 
product by p gives 

2 

pX-pQ-^—Q\X)=Q (22) 
u{p) 

Putting 

Q°(X) = e-^QPe** (23) 
the previous equation reduces to 

_ i ± + J* ?_Q0 = o (24) 

dp 2uj 2 (p) uj{p) 



Observe how, when acting on the detected subspace, Eq. (|21|) reduces ef- 
fectively to only the one-dimensional equation (Pl|). One would be tempted 
to solve it with the ordering chosen in (|12"D , with eigenfunctions similar to 



(Hj). This choice would not do, as the norm of these states would be badly 
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divergent in three dimensional space. What we need are eigenstates with 
higher negative powers of k than in ([14]). This can be achieved by choosing 
a different ordering for the operator. Tentatively we put 

^-^M-% + ^) p ^ ] (25) 

with this choice we get for the eigenfunction of fl23|) with eigenvalue T the 
expression 

tfrW(jfe) = J_^(k)T~kX) (26) 

where we have chosen some arbitrary fixed X. We now choose n such that 
the scalar product be well behaved 



2u(k) 

poo u2{l-n) 

= (27T)- 1 / dk^—-e^ T - T ^ (27) 
Jo uj(k) 

We see that the eigenfunctions are not orthogonal. We will address this 
problem in the next section. Now, we focus on the last integral, which 
strongly suggest the choice n — 1/2. In the general case of d space dimensions 
we would chose n = (d — 2)/2, to make the measure of the integral equal to 
duj. Finally, in our case we have: 

d ip 



¥ T X \k) = (2n)- 1 k- 1 / 2 e *C<"(*> r -**) (28) 

(4*),4f)) = (27T)- 1 r du,<Mk)<r-T) 



If there is any doubt left in that the right choice is n = 1/2, one can check that 
this value gives the unique ordering that makes the operator Q° Hermitian, 

(0,QV) = (QV,^)- 
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5 Orthonormalization and completeness 



The eigenfunctions of (28) are not yet orthogonal. However the above scalar 
product is an appropriate expression for the Marolf 's orthogonalization reci- 
pe [||. It is based in the physical observation that for vanishing momentum 
the particle either never reaches the detector, or sits in it forever. To deal 
with this situation, Marolf proposed a regularization prescription for the 
time-of-arrival operator that "avoids" zero momentum particles. The proce- 
dure to follow is less obvious here than in the 1-D non-relativistic case, due 
to the more complex structure of the operator. We first present the appro- 
priate prescription for arbitrary n, coming back to n — 1/2 at the end of the 
calculation, to show that only with this value the procedure gives orthog- 
onal eigenfunctions in three space dimensions. First, we rewrite Q° in the 
momentum representation as 

1 d b n+1 / 2 

Q° = —iu(k) t—±L— (29) 

which we regularize as follows 

Q° = -Mk)^)^ 2 ±k^^), (30) 
and where / is the same as in [1] 

M = {UkZt<l (31) 

The eigenfunctions ^j^\k) corresponding to this operator are of the form: 

{ x),^ 1 e^ T -^ _ r k dk' 



VP>(k) = ==, Z(k) = / , r , r . (32) 

1 } ^k^y/Jik) I* u(k>)f(k>y 1 ; 

and the orthogonality condition reads 

(ib {X) il/ X) ) - (2tt)- 2 f ^ k 1 r lZ(k)(T '- T) > (W) 

* } - {2n} J 2uj{k)f{k)k^ 6 {66} 

For the case n = 1/2 one gets 

(#), )) = (27T)- 1 f ZmaX dZe^ T '- T ^ = 5{T - T'), (34) 
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as the coordinate Z goes from —oo to as A; goes from to e, and from to 
oo as A; goes from e to oo. Z and T form a pair of "conjugate" variables in the 
subspaces H^q. This can be seen from ([53]) and the associated completeness 
relation 

£J dT*?\k)*P* (k') = —lj^ S{Z{k) - Z{k')) e*$-*# (35) 

The weird expression on the rhs is exactly what is needed to form a com- 
pleteness relation in the detected subspace. For any function G H^g 

I 2^) {Cy^^T^k')) = (36) 

as should be expected. In addition, using the expressions (|30| ) for Q° and 



(32) for Z, the following commutation rule is derived 

[Q°, Z] = -1 (37) 

The spectral support of both Q° and Z is the whole real line, so that no dif- 
ficulties arise from the Stone - Von Neumann theorem with (|37| ) as would be 
the case were it to involve uo instead of Z. Finally, a comment on the relation 
between the time and the position operators is in order: The eigenstates of Q 
with eigenvalue X (||) belong to the detected subspace H^q. However, it is 
not possible to determine simultaneously both the position (or the momen- 
tum) and the time-of-arrival due to the fact that the corresponding operators 
do not commute. 



6 Interpretation 

The results obtained so far indicate that the operator formalism associated 
to the time-of-arrival at a point works to fit the quantum mechanical rules. 
Accordingly, one can interpret it in a novel but standard way as was done 
on physical grounds in ref |TJ for one space dimension. Here, we will show 
that the formalism provides the tools with which to build the quantum me- 
chanical interpretation to be given to the time-of-arrival operator. In other 
words, that it provides the mathematical framework sufficient to define the 
time-of-arrival properties of the particle and associate to them definite proba- 
bilities. For definiteness, we assume that we are analyzing the time-of-arrival 
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at the point X. First, we split the Hilbert space Ti of states into never de- 
tected Hnd and detected subspaces Hd', obviously H = Hd © 7~Cnd- Also, 
from the discussion in Section 4, we know that H D = H {x) . This will be 
the Hilbert space appropriate to the analysis. In we have defined the 
(regularized) Hermitian operator Q°(X), whose spectrum is T e 1Z, the set 
of observable times-of- arrival at the point X. Having solved the eigenvalue 
problem for Q°(X), we obtained a complete and orthogonal set of eigenfunc- 
tions ip T (k) =< k\T,X > in the momentum representation. From them, 
we can define the set of elementary projectors {n^^T e 1Z} where 

= \T,X >< T,X\ (38) 
They generate a boolean algebra B with the properties 

nP f = n<r°, n^n^ = 8(t - t^u^ (39) 

To each elementary projector there corresponds an event (n^ <-> arrival at 
time T). Given any two projectors 11,11' E B the meet (and) and join (or) 
operations are defined as usual by 

n a n' = nn', nvn' = n + n'-nn' (40) 

where the notation corresponding to a finite dimensional Boole algebra has 
been displayed for simplicity. Statements will in general be of the form 
(Q°(X),Ti < T < T 2 ), i.e. the particle arrives at X in the interval (T 1: T 2 ). 
Associated to them there will be projectors built by the joining of elementary 
projectors of the algebra 

n (x) (T\, T 2 ) = r dT Tl^ (41) 

with matrix elements 

< T, X|n (x) (Ti, T 2 ) \T', X >= 5{T - T')6(T 2 - T)9(T - T x ) (42) 

Finally, the algebra has to provide a decomposition of the identity suitable 
for the analysis of the properties of the observable under discussion, i.e. 

n (x) = f + °° dTYl { T ] = 1 (43) 
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which is valid in due to with the obvious meaning that an arbitrary 
state oin (x) will not escape from detection. When acting on states belonging 
to Hilbert spaces larger than the value of will be smaller than one. 

The complement of the statement U ( - x \Ti,T 2 ), i.e. the particle arrives 
at X at a time outside the interval (Ti,T 2 ) will be given by the projector 
Y[( x ) — H( x }(Ti,T 2 ). In the case that the state of the particle belongs to 
the complement gives simply 1 — IT^ (Ti,T 2 ). The statement that 
there are states that escape from detection, absolutely when their projection 
on the detected subspace vanishes, or partially when they do not belong to 
liS x ^ but have a finite projection on it, is given by the projector 1 — ~^ x ^ . 
Finally, joining this last to the complement, gives the negative statement 
1 — n( x )(Ti, T 2 ), i.e. the particle does not arrive at X in the interval (Ti, T 2 ). 
The fact that the negation and the complement may differ is a consequence of 
the incomplete character of the spectral decomposition of the time-of-arrival 
operator (H^ x ' < 1). This could be avoided by working inside HS X > only, but 
this is too small to be of practical interest, consisting only of spherical waves 
about X. 

We can now assign probabilities to the statements represented by the 
projectors of the algebra B. Given an arbitrary normalized state $ of the 
physical Hilbert space, the probability (in time) of arriving during the interval 
(T 1; T 2 ) at the position X, P^\&) is given by 

p£ T #) = Q dT I < T '^l $ > I 2 ( 44 ) 

An arbitrary state $ does not need to be in Tin, but in general will have a 
finite projection on it. Accordingly, we can define the probability of being 
ever detected at X by 

pW($) = / dT\<T,X\$>\ 2 (45) 

This will be equal to one for normalized states in Tip, as can be obtained from 
(HD- For states not in Hd this describes the case of states that classically 
would never be detected at the position X, but quantum mechanically have 
a -less than one, but finite- probability for (ever) being detected at that 
point. Consider for example the ideal situation in which we place a detector 
along the ox axis at X = (x, 0,0), and prepare at t = a gaussian wave 
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packet centered at the origin, with mean momentum slightly off the ox axis 
< k >= (ko sin 9, 0, ko cos9). We consider the uncertainties in position and 
momentum to be such that wave packet and detector are well separated at 
t — 0, and the cone of flight of the particle (89 ~ Ak/k) misses the detector. 
Even in this case, there will be a small probability for the particle being ever 
detected at X; it is given by p( x )($). The probability of being detected 
during the interval (Ti,T 2 ) will be given by ^ > (V 1 t 2 )(^)' wn il e the average 
value of the time-of-arrival operator will be 

n,-x f + ^dTT\ < T,X\$ > I 2 , x 

< Q\X) >= ±£ 1 ^ L 46 

f + ™dT\< T,X\$ > I 2 

This is a conditional average value, i.e. it makes sense only in the case when 
the particle is ever detected. Speaking about the value of the time-of-arrival 
in the other case is a logical contradiction, undefined mathematically, as in 
this case < T,X\$ >= 0. 

The question of the time-of-arrival still deserves further clarification in 
quantum mechanics. We have outlined the mathematical framework whose 
existence allows for the asignment of probabilities to its different statements 
and for the use of logic to make inferences. In doing this, we are implic- 
itly considering the existence of measurement devices (detectors in this case) 
which will function almost ideally, without introducing serious disturbances 
in the experimental results, so that the logical outcomes can be compared 
straightforwardly with the actual results. The existence of such detectors 
goes beyond the scope of the present work, which only deals with the formal- 
ism and its interpretation. This is a question common to this (distributions 
in time), and the usual (distributions in space) formulations of Quantum 
Mechanics, and we can think that what is applicable there is also applicable 
here. Other serious issue, of actual interest for its practical implications, is 
the inclusion of interactions in the formalism. For instance, How will the 
gravity field of the Earth modify the distribution of times-of-arrival as mea- 
sured on the laboratory? This is of interest as there are experiments based 
on the production of a time-of-flight spectrum against the force of gravity. 
Another question is that of the time-of-arrival at a detector of a particle 
after traversing a barrier by quantum tunnelling. There is no classical ana- 
log to this situation. Therefore the method presented here will be useless 
to address this problem, which calls for a completely quantum mechanical 
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approach. There is a long list of pending questions worth of further research. 
Here, we turn to one of the motivations of this work: using the relativistic 
particle as a guideline to learn about time in quantum gravity. In principle, 
it would be plausible to think of the space part of the metric as playing a 
role similar to that of the detector position. Then, constraints restricting the 
detected Hilbert space as in (fl5|) are likely to appear. Were this the case, the 
comparison would be among different possible initial states (of the Universe 
(?)), and the subject of comparison the time employed by these states to -or 



the probability of- "evolve" |TT[] to a definite space metric. All this is highly 



speculative and object of further research. First of all, it is not even clear 
the mere existence of a suitable classical scheme from which to derive a time 
operator in the general case. 
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